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1 .  Introduction 


The  problem  of  active  feedback  stabilizing  flexible  structures  is  of 
recent  interest  especially  for  deployment  in  space  [1] .  A  significant  feature 
of  this  application  is  the  need  for  "robust"  controllers  -  whose  design  does 
not  require  precise  knowledge  of  system  parameters  —  see  [2] .  Although  a 
class  of  such  linear  stabilizing  control  laws  is  given  in  [2],  it  would 
appear  that  to  generate  the  control  effort  necessary,  the  actuators  (reaction 
jets,  for  example),  would  need  to  be  nonlinear:  linear  for  small  amplitudes 
but  nonlinear  (saturating- type)  for  large  amplitudes. 

Another  feature  of  the  space  application  is  the  uncertainty  in  the  damping 
parameters  (even  including  the  model) .  The  strong  stabilizability  using 
linear  controllers  requires  controllability  and  this  was  established  in  [2] 
neglecting  damping;  but  the  proof  of  controllability  when  damping  is  present 
is  unclear.  In  this  paper  we  show  that  for  a  reasonable  damping  model  the 
robust  linear  controller  still  yields  strong  stability  under  a  natural 
extension  of  a  condition  involving  the  undamped  modes  used  in  [2].  Under  the 
asumption  that  there  is  nonzero  natural  damping,  however  small,  we  prove  that 
we  have  strong  stability  for  a  class  of  nonlinear  feedback  control  laws  which 
are  linear  for  small  amplitudes  and  can  saturate  for  large  amplitudes. 

We  only  treat  the  abstract  wave-equation  version  of  the  problem:  the 
reduction  of  the  physical  model  (biharmonic  beam  equations  with  delta-function 
controls)  to  the  abstract  form  is  given  in  [3].  We  only  note  here  one  feature 
of  the  problem:  the  differential  operators  do  not  involve  any  "boundary 
conditions":  and  although  the  control  is  exerted  on  the  boundary,  the  formula¬ 


tion  does  not  fall  into  the  class  of  boundary-control  problems  treated  by 


Lasiecka  and  Triggiani  [4]  so  that  in  particular  the  results  therein  are  not 
directly  applicable. 

Finally  we  note  that  our  results  may  be  regarded  as  an  extension  of 
the  Benchimol-Slemrod  result  [5]  (see  Levan  [  6]  for  a  recent  treatment)  to 
a  class  of  nonlinear  controllers  but  without  invoking  controllability  —  albeit 
in  a  particular  case.  Our  proof,  although  totally  elementary,  relies  heavily 
on  the  eigenfunction  decomposition  of  the  generator.  The  problem  statement 
and  the  main  results  are  in  Section  2. 

2.  General  Results 

Let  H  denote  a  separable  Hilbert  space  and  let  us  consider  the  follow¬ 
ing  canonical  abstract  differential  equation  characterizing  the  response 
x(t)  of  a  flexible  structure  to  the  applied  input  u(t): 

Mx(t)  +  Dx(t)  +  Ax(t)  +  Bu(t)  =  0  (2.1) 

where  M  is  a  self-adjoint  positive  definite  (zero  is  in  the  resolvent  spec¬ 
trum)  linear  bounded  operator  mapping  H  into  H; 

A  is  a  self-adjoint  nonnegative  definite  closed  linear  operator  with 
dense  domain  and  with  compact  resolvent;  we  shall  (for  simplicity)  also 
assume  that  zero  is  in  the  resolvent  set. 

D  is  a  self-adjoint  nonnegative  definite  closed  linear  operator  whose 
domain  includes  that  of  ^A  ;  and 

B  is  a  finite-dimensional  linear  operator  mapping  r"'  into  H. 
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The  next  step  is  to  introduce  the  "energy  norm"  space,  which  we  shall 
denote  by  H  .  On  the  product  space 


P(/a  )  X  H 


we  can  introduce  the  "energy"  inner  product: 


(2.2) 


(Y,Z]g  =  [/a  y^^,  /a  z^]  +  [My^,  z^] 


(2.3) 


where 


Y  = 


Z  = 


(2.4) 


and  complete  the  space  in  this  inner  product  to  yield  To  avoid  confu¬ 
sion  we  shall  use  a  subscript  E  to  denote  the  inner  product  in  Let 

A  denote  the  operator 


A  = 


-m'^a  -m“^d 


(2.5) 


with  domain: 


- 

y. 

y,  e  P(A)  1 

V(A}  = 

Y 

L 

Y  = 

j. 

^2 

'  ^  V[/h) 

(2.6) 

We  shall  now  show  that  A  is  closed. 


^1  ,n 

=  Y  e  i:>(A)  ; 

y  Y  = 

^1 

^2,n 

n 

n 

^2 

where 


=  y2,n 


(2.7) 


Mz  =  -Ay,  -  Dy_ 
2,n  l/H  ■^2,n 


(2.8) 


and 


l"n- 


II '/a  (z  -  z  )  I  +  [(M(z  -  2  ),  (z 

ifH  A  l^n  n  2^11 

Since  (i^ )  has  a  bounded  inverse,  it  follows  that 


converges.  But  by  (2.7) 


z  =  (/a  )  ^  (/a  z  ) 

l,n  l,n 


y2,n  =  "l,n 


so  that 


{y-  }  and  {/K  y  } 

z  ,n  z  ,n 


are  Cauchy.  Now  since  domain  of  D  includes  that  of  /a 


D(»^  )' 


-1 


is  linear  bounded,  and  hence 


Dy  =  D(/a  )  ^ (/a  y  ) 
2  ,n  2  ,n 


converges.  Since  the  left  side  of  (2.8)  is  Cauchy,  this  implies  that 


1  ,n 


is  Cauchy,  and  of  course  {y,  }  is  Cauchy.  Hence  we  see  that  A  is  closed. 

1  ,n 


Since 


with  same  domain  as  A,  and  for  Y  in  P(A); 


■2lDy^,  y^] 


where 


[Ay,  Y]^  +  [A*y,  Y]g 


Y  = 


we  see  that  A  and  A*  are  dissipative,  and  hence  [ 7  ]  A  generates  a 
strongly  continuous  contraction  semigroup  which  we  shall  denote  S(t), 
t  ^  0.  Let  Aq  denote  the  "undamped"  generator: 


-M“^A  0 


(2.9) 


with  domain 


Y  = 


fi  P(A) 


Let 


V 


0  0 
0  M~^D 


We  note  that  A  >  0  belongs  to  the  resolvent  of  A^  and  that  R(A,Aq) 
is  compact.  From  the  easily  verified  resolvent  equation: 


R(A,  A)  =  R(A,  A^)  [I  +  PR(A,A)]  (2.10) 

it  follows  that  R(A,A)  is  compact,  since  PR(A,A)  is  bounded. 

Rewriting  (2.1)  as 

Y(t)  =  AY(t)  +  8u(t)  (2.11) 


I 
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where 


Y(t) 


x(t) 

x(t) 


Bu 


0 


we  see  that  (2.11)  has  the  "mild"  solution: 

t 

y(t)  =  S(t)  Y(0)  +  /  s(t-o)  8u{a)  do  (2.12) 

0 

Let  Sq(-)  denote  the  semigroup  generated  by  A^.  Let  denote 

the  eigenvectors  of 

M-^A 

orthonormalized  so  that 

2 

and  such  that  are  monotone  increasing.  Note  that  these  are  the  undamped 

or  "natural"  modes. 

We  shall  say  that  a  (time-invariant)  "feedback"  control 


u(t)  =  K(Y(t)) 


(2.13) 


where  K(-)  maps  into  R  ,  "stabilizes"  the  system  (2.11)  if 


Y(t)  =  S(t)  Y(0)  +  /  S(t-o)  BK(Y(a))  do 

0 

has  a  unique  strongly  continuous  solution  such  that  it  is  globally  stable. 
That  is  to  say 
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I 


|Y(t)ll^  -  0 


as  t  -► 


for  every  initial  y(0). 

We  begin  with  linear  controllers. 


THEOREM  2.1.  Suppose 


(D  +  BB*)(t)j^]  ^  0 


(2.14) 


for  every  k.  Then  the  linear  feedback  control 


u(t)  =  -p5*Y{t) 


(2.15) 


where  P  is  positive  definite  and  has  a  bounded  inverse,  yields  global 
asymptotic  stability. 

Conversely,  (2.14)  is  necessary  if  these  controls  are  to  yield  global 
stability. 

Proof.  To  begin  with  we  assume  that  the  strong  version  of  (2.11)  is 


satisfied: 


Mx(t)  +  Dx(t)  +  Ax(t)  +  Bu(t)  =  0  . 


(2.16) 


E(t)  =  i  iy(t);i^ 


where 


=  ■|([/Ax(t),  /Ax(t)]  +  (Mx(t),x(t)]) 


—  E(t)  =  -[Dx{t),  x(t)l  -  [x(t),  Bu(t) 


=  -[ (D  +  A(t)BB*)x(t) ,  x(t) 


[u(t) ,  B*x(t) 

II  _  ^  II  2 


B*x(t) i 


(2.17) 


(2.18) 


(2.19) 


(2.20) 


if  l|B*x(t)||  =  0 


.V'.t  >  •  >  .  .. 


Using  this  expansion  in  (2.16)  we  have 


or. 


aj^(t)  +  Uj^aj^(t)  =  -[Dx(t), 


[B*x(t)  ,  B*((i|^] 


(2.26) 


(2.27) 


where 


f^(t)  =  [Dx(t),4ii^]  ~  (t)  +  [B*x(t)  ,  B*((i 


aj^(t)  [B*(J)j^,  B*(})j^] 


Suppose  now  that  for  some  k: 


D\  =  0  '• 


B*())  =  0 

iv 


Then 


x(t)  =  aj^(t)iiij^ 


is  a  solution  of  (2.17),  provided  only  that 

2 

aj^(t)  +  Wj^aj^(t)  =  0 

and  hence  ^  does  not  go  to  zero  as  t  -*■  this  takes  care  of  the 

necessity  condition.  Assume  then  that 


20]^  -  1%. 

is  nonzero.  Then  (2.27)  can  be 
for  simplicity) : 


"solved"  to  yield  (assuming  small  damping 


ak(t)  =  e 


cos  X, 


+  b,  e 
k 


-Okt  sin  X^t 


t 

/  W(t-s)  f  (s)  ds 
0  ^ 


(2.28) 


where 


e  sin  A,  t 
W(t)  =  - r - ^ 

\ 


(-0,  ±  iX, )  are  the  roots  of 
k  k 


2  .  2  „ 

s  +  2o,  s  +  0)  =  0 

k  k 

Now  because; 

|iDx{t),  S  [Dx(t),  x(t)]  [D())^,  4.j^] 


it  follows  from  (2.23)  that 


/  I  [Dx  (t)  ,  (|),  ]  I  ^  dt  <  “  . 

0  ^ 


Hence  a  little  analysis  shows  that 


t 

/  W(t-s)  [Dx(s),  (}>  ]  ds  0  as  t  -*■  ” 

0  ^ 


By  a  similar  reasoning,  also 
t 

/  W(t-s'  [B*x(s),  B*'(})  ]  ds  -»■  0 

0 


as  t 


Hence  it  follows  that 


E(t)  =  a,  (t)  -  2o  /  W(t-s)  a  (s)  ds 

k  k  _  k 


where 


e(t)  -  0 


as  t 


By  integration  by  parts  we  obtain  that 


=  +  3^(0)  W(t)  +  2o^  /  W(t-o)  3)^(0)  do 


e  (t) 


Hence 


t 

a  (t)  +  20  /  W(t-a)  a  (0)  da  =  6(t)  (2.30) 

K  R  Q  K 

where 

6(t)  -*■  0  as  t  -►  ®  . 

Solving  (2.30)  (using  elementary  Laplace  transform  techniques)  we  obtain 
that 

a,  (t)  ->•  0  as  t  -►  *  . 

k 

Differentiating  (2.28)  and  proceeding  in  a  similar  way  we  can  also  show  that 


Let 


Then 


for  each  )<,  which  is  enough  to  imply  that  Y(t)  converges  wea)cly  to  zero. 
Next  we  note  that 

Y(t)  =  S„(t)  Y(0) 

O 

where  S  (t)  is  the  (strongly  continuous  contraction)  semigroup  generated 
B 
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which  has  a  compact  resolvent.  Hence  weak  stability  implies  strong 
stability  (  7  ] : 


|Y(t) 


|Sg(t)  Y(0)ij. 


0  . 


(2.31) 


V 


In  particular  we  have  (2.23)  that 


|y(0) 


/  [Dx  (t)  ,  X  (t)  ]  dt 
0 


/  ||B*x(t)  II  ^  dt 
0 


estimate : 


Y{0) 


=  /  [DY(t),  Y(t)]  dt  +  /  ||B*Y(t) 


dt 


(2.32) 


Finally  given  any  arbitrary  initial  condition  Y(0),  we  can  find  an 

approximating  sequence  {Y^(0)}  in  the  domain  of  A  such  that  (2.31)  holds 

for  the  corresponding  solution  Y  (t) ;  and  the  result  follows  from  the 

n 


|iY(t)iL  ^  l|S„(t)  Y  (0)|L  +  IIy  (0)  -  Y(0) 


RE.'IARK  1.  We  note  that  if  D  =  0,  our  condition  (2.14)  is  equivalent  to 
requiring  (A ~ B)  controllability  (see  (2]).  Hence  we  would  have  strong 
stability  by  a  more  general  argument  due  to  Benchimol  iS).  If  B  =  0,  then 
condition  (2.14)  implies  that  Di}>  /  0  for  any  k,  and  we  are  then  proving 

Iv 

strong  stability  for  the  semigroup  S(t). 

Nonlinear  Controllers 

Let  us  now  go  on  to  consider  a  class  of  nonlinear  controllers.  Thus 
we  shall  consider  where  K(')  in  (2.13)  is  given  by 
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K(Y)  = 

where 

Y  = 

and  f ( • ) 

maps  R  into  R  , 

i) 

[f (u) ,  u]  >  0  for 

ii) 

l|f(u)ll  S  Xiui 

iii) 

f  ( • )  is  Lipschitz. 

A  typical 

exam.ple  of  f  ( • )  is 

f (B^y^) 


and  satisfies  the  following  conditions: 

M  0 


f(u)  =  V  ; 


u  =  { u .  }  ; 


V  =  {v. } 


V .  =  Y .  tan  ^u.u.  ;  Y.»U.  >  0 

1  1  11  11 


THEOREM  2.2.  Suppose 


(2.14a) 


for  every  k.  Then  the  feedback  control 


u{t)  =  -B*  f (B*x(t) ) 


(2.33) 


where 


Y  (t) 


x(t) 

x(t) 


yields  asymptotic  global  stal ility. 

Proof.  Under  condition  (2.14:',  the  semigroup  S(')  generated  by  A 


in  H  is  strongly  stable.  By  virtue  of  the  Lipschitz  conditions  on  f(*)/ 
existence  and  uniqueness  of  solution  for  each  t  >  0  is  immediate,  and  the 
solution  is  given  by: 

t 

Y(t)  =  S(t)  Y(0)  +  /  S(t-O)  B  f (B*x(o) )  do,  0  <  t  .  (2.34) 

0 

Next  we  shall  show  that 


/  ||f  (B*x(a))  do  <  »  .  (2.35) 

0 


For  this  purpose  let  us  assume  first  that  the  initial  condition  Y{0)  is  in 
the  domain  of  the  generator  A,  so  that  we  have 


Mx(t)  +  Dx(t)  +  Bf(B*x(t))  +  Ax(t)  =  0  . 

Defining  the  energy  again  as 


we  have  that 


where 


Since 


we  have  that: 


E(t) 


E(t)  =  I  ||Y(t)jlg 
-(Dx{t) ,  x(t) ]  -  [b(t) ,  f (b(t) ) ] 


b(t)  =  B*x(t)  . 


(b(t),  f(b(t))]  S  0  , 


/  [Dx (t) ,  X (t) ]  dt  <  “ 

0 


.1  .t  ^  *,»  *  *  '**  'k*  *  • '  I  I  *•»  * '  * » *  . 
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I 


and  since  we  can  write 


it  follows  that 


But  this  implies  that 


and  hence  also 


(/d  (/Dx(t))  =  x(t) 


/  ||x(t)  dt  <  “  . 

0 


/  i|B*x(t) ir  dt  <  “ 
0 


J  if(B*x(t))i  dt  <  » 

0 


by  virtue  of  our  assumptions  on  f(*). 

We  now  proceed  as  in  the  proof  of  the  preceding  theorem.  Writing 


we  have ; 


♦ki 


*  Vk"’  *  2Vk‘’=’  '  -'k*'’ 


(2.36) 


where  a,  (t)  is  as  before,  and 
k 


fj^(t)  =  [Dx(t),  -  aj^(t)  +  [f(B*x(t)), 


/  [Dx(t),  <{)|^]  dt  +  J  [f(B*x(t)),  dt  < 


as  before,  we  obtain 

t 

e(t)  =  w(t)  +  20^^  /  w(t-o)aj^(a)  do 


i'w  I*  V  .  :* 


S' 


where  e(t)  -*■  0  as  t  “  and  the  rest  of  the  arguments  follow.  Hence 

aj^(t)  0  as  t  “ 

aj^(t)  -*•0  as  t  -►  ® 

for  each  k.  Hence  it  follows  that 

Ax(t)  0 

Ax(t)  0  . 

Since  A  ^  is  compact,  it  follows  that  x(t)  converges  strongly  to  zero 
(in  H) .  Similarly 

Ax(t)  0 

implies  that  x(t)  converges  strongy  to  zero  (in  H) .  Hence  Y(t)  converges 
weakly  (in  H^)  to  zero.  Now 


implies  that 


Hence 


(Dx  (t)  »  =  [x  (t)  , 

Dx(t)  0  . 


AY(t)  ^  0  . 


since  A  has  a  compact  resolvent  it  follows  that  Y(t)  converges  strongly 
(In  Hj.)  to  zero. 

Next  we  relax  our  assumption  regarding  the  initial  condition  Y(0). 

We  begin  by  proving  Lipschitz  continuity  of  the  solution  with  respect  to 
the  initial  condition.  Thus  let  Y(t)  denote  the  solution  with  initial 


condition  Y(0)  and  let 


Y(t) 


M(t)  (Y(0)  )  . 


Then  let  Y,(0)  ,  YJO)  «  and 
1  Z  L 

M{t)(Y2(0))  -M(t)(Y^(0))  =  Sj^(t)(Y2(0)  -Y^{0)) 

t 

+  /  Sj^(t-o)  (Bf  (8*Y2{o)  )  -  Bf(B*Y^(a)))  do 

Let 

in(t)  =  !iY2(t)  -  Y^(t)|lg  . 

Then  in  view  of  our  assumptions  on  f('),  we  have 

t 

m(t)  S  m(0)  +  y/  into)  do  ;  Y  -  0 

0 


and  hence  by  the  usual  analysis 


m(t)  -  m(0) 

yielding  Lipschitz  continuity.  The  continuity  yields  in  turn 


’!Y(t)i^  =  2E(t)  =  I|y(0)||^  -  /  [Dx(5)  ,  x(5)  )  ds 


t 

/  [b(s),  f(b(s))]  ds  , 
0 


and  hence  we  obtain: 


/  iDx(t)li^  dt  +  /  [b(t),  f(b(t))]  dt  <  “ 

0  0 


Next  we  need  to  establish  (2.36).  But  this  follows  readily  from  the  fact 
(2.11)  holds  in  the  wea)c  sense.  Hence  we  obtain  that  Y(t)  converges 
strongly  to  zero  in  H  for  any  initial  Y(0)  in  H^. 

REMARK.  It  would  be  of  interest  to  establish  strong  stability  under  the 
wea)cer  condition  (2.14).  We  note  that  the  condition  (2.14)  is  again  obviously 
necessary . 
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